Weyl formulas for annular ray-splitting billiards 



o 
o 

(N 



O 

(N 



Q 

U 
d 



> 

m 
o 
in 
o 
m 
o 



X 



Yves DecaninQ and Antoine Folacc{j] 

&P-E, t/M7? CNRS 6134, Equipe Physique Semi-Classique (et) de la Matiere Condensee 
Universite de Corse, Faculte des Sciences, BP 52, 20250 Corte, France 
(Dated: February 8, 2008) 

We consider the distribution of eigenvalues for the wave equation in annular (electromagnetic or 
acoustic) ray-splitting billiards. These systems are interesting in that the derivation of the associated 
smoothed spectral counting function can be considered as a canonical problem. This is achieved by 
extending a formalism developed by Berry and Howls for ordinary (without ray-splitting) billiards. 
Our results are confirmed by numerical computations and permit us to infer a set of rules useful in 
order to obtain Weyl formulas for more general ray-splitting billiards. 
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I. INTRODUCTION 

Since the pioneering work of H. Weyl [l| dealing with 
the distribution of eigenvalues for the wave equation in 
a cavity with a perfectly reflecting boundary consider- 
able effort has been devoted to the construction of the 
smooth part of spectral counting functions (Weyl formu- 
las) in various fields of physics and mathematics. For an 
historical point of view on this problem, we refer to the 
seminal paper of Kac Q and to the monograph of Batles 
and Hilf For its importance in physics and for var- 
ious applications, we refer to the monographs of Batles 
and Hilf Q and of Brack and Bhaduri p[ as well as to 
references therein. 

In general, the determination of the smooth part of a 
spectral counting function is a complicated task and only 
the leading order terms can be obtained. By contrast, for 
two- and three-dimensional billiards with Dirichlet, Neu- 
mann or Robin boundary conditions, this problem can be 
considered as definitely solved (see, e.g., Refs. IEIS 000 
or the monographs cited above). For ray-splitting bil- 
liards introduced in the context of acoustic and quantum 
chaos by Couchman et al |1C| and extensively studied 
these last years (see, e.g., Refs. [H d d [j E E 
ITtI ITsl IT^) . it is rather natural to think that the same 
outcome could be obtained. Recently^some progress has 
been made in that direction jllL Il3l [14. flq but it seems 
we are very far from a general theory. With this aim in 
view, it is interesting to solve canonical problems, i.e., to 
consider simple examples of ray-splitting billiard prob- 
lems for which it is possible to perform exactly the cal- 
culations 0, E] or to carry on them as far as possible 
[l6| . The results then obtained are useful in order to infer 
Weyl formulas for more general ray-splitting billiards. 

With this in mind, we are concerned, in this paper, 
with the distribution of eigenvalues for the scalar wave 
equation in a two-dimensional dielectric annular billiard. 
This billiard consists of an outer circle with radius R 



and an inner circle with radius r. The index of refraction 
between the two circles (region I) is fixed at 1 while the 
inner disk (region II) is characterized by the index of 
refraction N. At the interface between the two regions, 
we shall assume that the scalar field $ solution of the 
wave equation and its normal derivative d<fr/dn satisfy 
the boundary conditions 



and 



dn 



dn 



(1) 
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The cases a — 1 and a — 1/N 2 respectively correspond 
to the TM and TE polarizations in electromagnetism 
pfj|. On the outer circle, we shall assume that $ van- 
ishes (Dirichlet boundary condition). Such a condition 
is more artificial than physical. It can be partially real- 
ized for the TM polarization if the billiard is embedded 
in a perfect conductor. We assume it in order to simplify 
our calculations. Mutatis mutandis, we can also consider 
the case of a two-dimensional acoustic annular billiard. 
In that case, region I (resp. region II) is occupied by a 
perfect fluid with density p l (resp. p 11 ) while a — p 1 /p u 
|20j and can take any positive value. Finally, in order to 
be able to analytically perform the calculations, we shall 
assume that the two circles are concentric, but the final 
results given by Eqs. 1291311 and 0U| are equally valid for 
non-concentric circles. 

For such billiards, eigenvalues cannot be analytically 
obtained. They satisfy a transcendental equation involv- 
ing Bessel functions that can be solved only numerically. 
In spite of this, we are able to analytically derive the 
associated Weyl formulas from the corresponding Green 
functions. This is done by using an approach developed 
by Berry and Howls § and which generalizes a previous 
work by Stewartson and Waechter 7] . This approach has 
been considered by these authors for the circular billiard 
with the Dirichlet condition on its boundary. We extend 
it rather naturally to the more complicated case of an- 
nular ray-splitting billiards. We are then confronted by 
some tedious algebraic calculations which, fortunately, 
can be performed with the help of Mathematica [2l|. 

Our paper is organized as follows. In Section 2, we in- 
troduce our notations and we construct the Green func- 
tion for the annular ray-splitting billiard as well as the 
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associated regularized resolvent. In Section 3, by extend- 
ing the Berry-Howls approach, we obtain a set of Weyl 
formulas corresponding to various values of the parame- 
ters a and N. In Section 4, we briefly consider the same 
problem for the desymmetrized versions of the annular 
ray-splitting billiard. In Section 5, we numerically check 
the previous results and, in Section 6, we conclude our 
paper by inferring a set of rules useful for constructing 
Weyl formulas for more general ray-splitting billiards (see 
Eqs.EEHl). 



II. GREEN FUNCTION FOR THE ANNULAR 
RAY-SPLITTING BILLIARD AND 
REGULARIZED RESOLVENT 

From now on, we shall use the polar coordinate sys- 
tem (p, 9) with its origin O at the common center of the 
two circles which define the annular ray-splitting billiard. 
The eigenvalues ki for the wave equation in this billiard as 
well as the associated eigenfunctions $i are determined 
by solving the following problem: 

• (i) The ki and <E>i satisfy the Hclmholtz equation 



fl»*i(x) =fc?$<(x) x{p,0) 



where 



-A x 

-(1/JV 2 )A, 



for r < p < R, 
for < p < r, 



(2) 



(3) 



with the Laplacian A x given, in the polar coordi- 
nate system, by 



A x = 



dp 2 



ld_ 

pdp 



1_ cP_ 



• (ii) The satisfy the boundary conditions 
= r,0)=$?(p = r,0), 



dp 

for < 6 < 2tt 



(p = r,9) 



dp 



r,8), 



(4) 



(5a) 
(5b) 

(5c) 



Because a solution of © is expressible in terms of Bessel 
functions (2^ , it is easy to prove from Eq. [5] that the ki 
are the values of k which solve 



J m {kR) J m (kr) 
H^(kr) 



J m (Nkr) NaJ' m {Nkr) 



Hm\kR) 



= 



(6) 



for m £ Z. For a given to, they can be indexed by the in- 
teger n with n = 1,2,3,... and the corresponding eigen- 
functions are given by 



H^{k m ^R) H ^ m '"P>) e ^) 



where the A m ^ n are normalization constants. It should 
be noted that the eigenvalues corresponding to to ^ 
are twofold-degenerated because of the relation k m ,n = 
k- m , n which follows from the invariance of Eq.[S] under 
the change to — ► —to. 

The determination of the eigenvalues fcj permits us to 
construct the spectral counting function associated with 
the annular ray-splitting billiard. It is given by 

+ OG +OG 

M(k) = ^6(fc-fc l ) = £ e ( fc -*>»,») ( 8 ) 

i m—~oo n—1 

where denotes the Heaviside function. 

Eq. can be solved only numerically. As a conse- 
quence, the smooth part of the spectral counting func- 
tion Af(k) cannot be obtained directly from (JSJ. In order 
to accomplish this, it is more convenient to generalize 
the Berry-Howls approach p| (see also Stewartson and 
Waechter Q) by introducing the regularized resolvent 



g(s)= r r [G n (x,x,«)-Gj?(x,x,a)] pdpd9 + 
Jo Jq 

J* J [G I (x,x, s )-G|,(x,x, s )] pdpdO (9) 
where Go is the "free-space Green function" given by 
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G (x,x', S ) = — tf ( fl |x-x'|; 

j +oo 

r Yl / m(sp<)^ m (sp; 



2tt 



(10a) 



N 



G I I (x,x',s) = — K (Ns\x-x'\ 
Ztt 



N 2 

— ]T I m {N SP< )K m {N sp> )e"^ e - e \ 



m— — oo 



(10b) 

while G is the annular ray-splitting billiard Green func- 
tion solution of 

(i? x + s 2 )G(x,x',s) = tf(x-x') (11) 

and subject to the boundary conditions 

G\x,x',s) = 

for p or p = R and < 9, 9' < 2vr,(12a) 
G^x^s) = G n (x,x',s) 

for p or p = r and < 6, 6' < 2tt, (12b) 

dG l , dG u , 

(x, x , s) = a ~^-( x > x > s ) 

for por p' = r and < 0, 9' < 2n. (12c) 

Here, it should be noted that in order to construct g(s), 
we need the Green function G(x, x', s) only for x and x' 
lying in the same region of the billiard. 
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When |s| is large, g(s) has the asymptotic expansion 
(Weyl series) 

= E 5 ( 13 ) 
p=i 

and from the c p -coefficients we can obtain the large-fc 
asymptotic behavior for the spectral counting function 
in the form 

« r/T\ A 2 2ci fc (-1)" 
■A/ 0) = J" * + k + c 2 - - 7 TTT^ C2n+1. 

47T 7T 7T ' (n — -k) h 

(14) 

I 



Our theory does not provide the expression for the sur- 
face term A. We shall assume that A is the billiard total 
area weighted by the refraction index and given by 

A = tt(R 2 - r 2 ) + N 2 nr 2 . (15) 



In order to obtain the c p -coefficients, we must first 
solve the problem defined by Eqs. ^2 and El and then 
to perform the integrations in Eq. The solution of 
(|llf> and l|12|) can be constructed in terms of the modi- 
fied Bessel functions |22j and is given by 



G 



(x,xVs) = - £ 



I m (sR)D^\sr, N) 



Hfoo I m (sR)D m (sr, N) - K m (sR)D^\sr, N) 



K m (sR) 
I m (sR) 



D m (sr,N) 
D^isr^N) 



Im{sp<) - if m (sp<) 



i(6- 



N ] 

G n (x,x', S ) = — Irn(N SP< ) 



2 +°° 



2tt 



K, 



(AT \ I m (sR)E m (sr, N) - K m (sR)E£> (sr, N) 
I m (sR)D m (sr, N) — K m (sR)D)n > (sr, N) 



with p< = inf (p, p'), P> = sup(/j, p') and 

D m (sr, N) = K' m (sr) I m (Nsr) - Na I' m (Nsr) K m (sr), 
Dg\nr,N) = I' m (sr)I m (Nsr) - Na I' m {Nsr) I m (sr), 
E m {sr,N) = K' m {sr)K m (Nsr) - Na K' m (Nsr) K rn (sr), 
E£)(sr, N) = I' m {sr) K m (Nsr) - Na K' m {Nsr) I m (sr). 

This provides the expression for g(s): 



where 



with 



9(s) = —x E f" 1 ^ 



fm(s)=f!n ) (s)+f!n ) (s)+f!n ) (s) 



I m (sR) K m (sR) - I' m (sR) K' m (sR) 



2 R 2 



J_ 4(gg) An(fr, AQ - (gr, AQ 

S-R / m (si?) D m (sr, 2V) - K m (si?) (sr, JV) 



/i 2) M = ^ 2 



1 



7 m ( S r)X m ( S r)-4( S r)^(sr) 



iVa I m {sR)K' m (sr)-K m (sR)I' m (sr) 



■l' m (Nsr) 



^ I m ( s R) D m (sr, N) - K m (sR) D$(sr, N) 

1 / m (si?)X m (sr) -X m ( S i?)/ m (sr) 



1 



sr I m (sR) D m (sr, N) - K m (sR) D^\sr, N) 



I m (Nsr) 



(16a) 
i(e-e') 
(16b) 

(17a) 
(17b) 
(17c) 
(17d) 

(18) 

(19) 



(20a) 



(20b) 



4 



f(3)r<A - N 2 r 2 



1 + 



1 



I m {Nsr) K m (Nsr) - J^(JVar) K' m (Nar) 

I' m {Nsr) 



N 2 s 2 r 2 

I m (sR) K' m {sr) - K m {sR) I' m {sr) 



N sr / m ( SjR ) Dm ( sr) AT) _ K m {sR) (sr, N) 

2 ^ a I m {sR)K m (sr)-K m {sR)I m (sr) 



1 



N 2 s 2 r 2) sr Im ( sR ) Dm ( s ^ N )_ Km ( sR ) D $( sr)N ) 



Im(Nsr) 



(20c) 



By using the Poisson summation formula as well as the 
relation /_ m (s) = / m (s), we can write 



9(s) 



J2 / fm(s)e l2 ^ m dm. (21) 



/x— — oo 



It should be noted that (|21|l with / m (s) given by H19|) 
and H2Ufl provides an exact expression for g(s). 



III. FROM THE REGULARIZED RESOLVENT 
TO THE SMOOTHED SPECTRAL COUNTING 
FUNCTION 



and Uj, and v p are polynomials given in chapter 9 of 
Ref.il (Eqs. 9.3.10 and 9.3.14). It should be noted that, 
as for the circle billiard @, II] , the Weyl coefficients c p 
and therefore the smoothed spectral counting function 
come directly from the /i = 0-term in (|21l) . As noted by 
Berry and Howls ||| (see also |23, 13 )> the other terms 
are associated with the fluctuating part of the spectral 
counting function which could be obtained by carefully 
taking into account Stokes phenomenon for the asymp- 
totic expansions (1221) in the context of hyperasymptotics 



The \i = 0-term in i|?T|) now reduces to 



The large-|s| asymptotic behavior l|13|) of g(s) can now 
be found from l|21|) by replacing in Eqs. ED and [3U] the 
modified Bessel functions I m ,I' m , K m and K' m by their 
uniform asymptotic expansions (see Eqs. 9.7.8 - 9.7.10 of 
Ref. IH given by 



9{s) 



f* + oo 



f m (s) dm. 



(24) 



I m {z) 



2ir (to 2 + z 2 ) 4 



- exp(^ m (^)/2) 



-\-oo 

,2 , „2\i 



Up \t m {z)] 



mP 



(22a) 



7^(z) = 1 V exp(^ m (z)/2) 



z^/27^ 



><E 

p=0 

1 



mP 



(22b) 



2 (m 2 + z 2 ) 4 



r exp(-F m (s)/2) 



X ^-^f^', (22c) 

p=0 



mP 



KM 



tt (m 2 + z 2 ) i 



exp(-F m (z)/2) 



z 

x g(-i)%M^ (22d) 



Here 



p=0 



top 



F m (z) = 2(to 2 + z 2 )5 + 2m log 

(TO Z + Z A ) 2 



TO + (to 2 + Z 2 ) 2 

(23a) 

(23b) 



where 



/m(«)=4 1) (*)+/^(*) + /^(«) 



(2), 



mi 



with 



/i 3) (s) = 



L — ' mP 

0=1 



(25) 



Vm 2 + s 2 R 2 ^ 4 1 ) [ tm ( SjR )] 

~ 2 j ~ > l/oaj 



Vm 2 + ,s 2 r 2 +^ A^_[t m (sr), t m (Nsr)} 

a 2 



mP 



(26b) 

a 2 2-^1 



mP 



(26c) 



Here the functions Ajf 1 , A p 2 ^ and A p 3 ^ can be expressed 
in terms of the polynomials u p and v p and therefore can 
be explicitly obtained (see, below, Eq. I28|) . Then, by 
using Eqs. 1241261 we find the general expression 



■5 



1 f +oa 



VI + x 2 



xP 



dx 



r p- 



[ + °° VT+xS A(2) 
Jo ^ p 



c/^/l + x 2 ,x/^/N 2 + x 2 



dx 



-/ 

Jo 



+°° VN 2 + x 2 
xp 



4 3) \x/y/l + x?,x/y/N*+x? 



dx 



(27) 



In order to provide the terms in k and k° in the expression of the smoothed spectral counting function, we need 
the coefficients c\ and ci. They can be obtained, by performing the integrations in the previous equation, from the 
functions A^\ A^ and A^ which are explicitly given by 



4 1 

A? 

A 2 

A? 

4(3 



x/\Jl + x< 
xj\J\ + x< 



2 ( 1 + x 2)3/2' 



2(l + a; 2 ) 3: 



xj\J\ + x 2 , x/\/N 2 + x 2 
x/^l + x 2 , x/y/N 2 +x 2 
x/yjl + x 2 , x/^/N 2 +x 2 



2(1 + x 2 f/ 2 (l + x 2 ) [VTTx 2 + a VN 2 + x 2 \ 1 
[(1 - a)N 2 + (1 - aN 2 ) 2 x 2 + (a - l)(2a + l)iV 2 a; 4 + (a 2 - l)x 6 ] 
2{N 2 + x 2 )(l + x 2 f [VTTx^+aV^ + x 2 ] 2 
N 2 x aN 2 x 



+ 



131 ' 



(28a) 

(28b) 
(28c) 

(28d) 
(28e) 



xj\J\ + x 2 , x/\/N 2 + x 2 



2 (N 2 + x 2 f/ 2 (N 2 + x 2 ) [^/T+~x 2 + a ^WT 
[a{a - l)N 6 + (1 - aN 2 ) 2 N 2 x 2 + (1 - a) (a + 2)N 2 x 4 + (1 - a 2 )N 2 x 6 ] x 2 



2(1 + x 2 )(N 2 + x 2 ) 3 [VTTx 2 + a VWT: 



For the TM polarization (a = 1), we then find 



(28f) 



M(k) 
Af(k) 



„2/n>2 _2\ 



<„2/E>2 _2\ 



+ 



+ 



N 2 k 2 r 2 


kR 








+ 


4 




N 2 k 2 r 2 


kR 








+ 


4 





kr +- + ■■■ if N < 1, 
6 



2N (N 2 -l\ N 1 
~ l^V 2- ) ~ ~2~ 2 



(29a) 



fer+- + -- - if iV> 1. (29b) 
6 



For the TE polarization (a = 1/N 2 ), we then obtain 
Af(k) 



k 2 (R 2 - r 2 ) N 2 k 2 r 2 _ kR 
4 + 4 2 



2 N 2 2 /V 4 

X(1-7V 2 ) " 



7T 



1 — TV 4 - 1-iV 2 

7T l + iv 2 V 1 ' 2' y 



1 

T ~ 2 



kr +- + ■■■ if AT < 1, (30a) 
6 



AT(A:) = 



k 2 (R 2 -r 2 ) N 2 k 2 r 2 kR 
4 + 4 2 



2 Af /Ar 2 -1 
iv 



7T 1 + A" 2 



A' 2 



n 



JV 4 -1 i JV 2 -1\ AT 1 



ttN(1 + N 2 ) \ N 4 ' 2' N 2 J 2 2 



kr+ - + H if AT> 1. 

6 

(30b) 



6 



Finally, in the general case (a ^ 1), we obtain 
k 2 (R 2 -r 2 ) N 2 k 2 r 2 



M{k) = 



Af(k) 



kR 

4 4 ~T 

2 a(l ~N 2 ) , , x 2 a 2 N 2 - 1 
- 2 ' K{l-N 2 ) + - n 

7T Or — 1 7T a(Q! z — 1) 



^,1-JV 2 



k 2 {R 2 -r 2 ) N 2 k 



2 q(l-JV 2 ) 
tt iV(a 2 - 1) 



A" 



kR 

4 IT 

'TV 2 - r 



JV 2 



2 a a 2 TV 2 - 1 / 7T N 2 
^—r -s II 1 - or 

7T JV(a 2 - 1) 



JV 



JV 2 



fer + - 
6 



TV 1 
T ~ 2 



fcr 



if JV < 1, 



(31a) 

if JV > 1. 
(31b) 



In order to perform the integrations leading to Eqs. I29t 
1311 it has been necessary to separate the particular case 
a = 1 (TM polarization) from the general one a =/= 1. 
It should be noted that the case a = l/N 2 (TE polar- 
ization) is included in the general case a ^ 1 and can 
be recovered from Eq. 1311 Moreover, it is important to 
keep in mind that, as far as the elliptic integrals E, K 
and II are concerned, we adhere with the definitions and 
conventions of Refj22j which are in agreement with those 
of Mathematica |2lj but differ to those of Ref. f27L 

The various terms appearing in Eas. I29l31l have their 
usual physical interpretations: The first and the second 
terms (in fc 2 ) yield the area contributions, the third and 
the fourth ones (in k) yield the perimeter contributions 
while the fifth term (in fc°) yields the curvature contri- 
butions. In particular, it should be noted that the fourth 
term, in all these equations, provides the perimeter cor- 
rection associated with the circular ray-splitting bound- 
ary at p = r and that it is this term which contains 
the elliptic integrals. Moreover, it seems to us necessary 
to point out that the fifth term is associated with the 
curvature of the circular Dirichlet boundary at p = R. 
In other words, and this is rather surprising, the ray- 
splitting boundary at p = r does not provide any correc- 
tion to the curvature contributions. 



IV. THE DESYMMETRIZED ANNULAR 
RAY-SPLITTING BILLIARD 

In Section 2, we pointed out the twofold degeneracy 
of the eigenvalues k m ^ n with m =/= 0. It is possible to 
work with the non-degenerated spectra by separating the 
eigenfunctions of the annular ray-splitting billiard in two 
different sets: In the first set, we consider the even eigen- 
functions (even in the change x(p, 9) — > cx(p, —9)) given 
by 



&+l(p,6) = AtiUj m (k m , nP ) 

^^H^p))^)^ 



with in 6 N while, in the second set, we consider the odd 
eigenfunctions (odd in the change x(p, 0) — > crx(p, — #)) 
given by 

^n{P' e ) = A ml( J m(km,np) 



with to G N*. Here and in the following the superscripts 
(+) and (— ) refer respectively to positive and negative 
parities and the Am} n and the Am} n are normalization 
constants. The spectral counting functions (fc) and 
Af 1 --^ (fc) associated with these two sets are then given by 



AA+>(fe) = J2 J2 e ( fc - fc ™.")' ( 34a ) 

m— n— 1 

+ OC +OC 

^ H w = EE e ( fc -M- (34b) 



m— 1 n— 1 



The smoothed spectral counting functions Af^(k) 
and Af(~\k) respectively associated with 7V' + '(fc) and 
J\f(~^(k) can now be obtained by using, mutatis mutan- 
dis, the theoretical framework developed in the two previ- 
ous sections: Af^* 1 (fc) can be constructed from the even 
part g( + \s) of the regularized resolvent g(s) given by 
@ while A/"(-)(fc) can be constructed from its odd part 
'(s). The functions g( + \s) and g^~\s) are given by 



3 (±) (s) = \f / 2X [G n (x,x, S )±G n (x,ax, S ) 
1 Jo Jo 

- G^(x,x,s)] pdpd9 
+ if P [G I (x,x, S )±G I (x,ax, S ) 

- G^(x,x,s)] pdpd9 (35) 

and they satisfy 

g(s)=g(+\s)+g^(s). (36) 
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By performing the integrations in (|35[1 . we obtain 

g {± Hs) = lg(s)±g COII (s) (37) 

with 

.9 corr (s) = ~\fo(s) 

N 2 r 2 

[I (Nsr)K (Nsr) - I^(Nsr)K' (Nsr)] 

[I (sR)K (sR) - I^(sR)K' a { s R)} 

-- [I (sr)K (sr) - I' ( S r)K' (sr)} (38) 

and, by using the as ymp totic expansions given by 
Eqs. 9.7.1 - 9.7.6 of Ref. HI we can write 

H.) = H^ -i + O (1). (39) 

4s 8S^ |s|->+oo \S J / 

We then immediately obtain 

^(k) = lm± iN ~l )r + R kTU- (40) 

with J\f(k) which is given by any of Ecis. I29I3"T1 according 
to the physical problem considered. 

Now, we would like to provide a physical interpreta- 
tion of the results obtained above. We first note that the 
twofold degeneracy of the eigenvalues k m _ n with m =/= 
is directly linked to the invariance of the annular billiard 
under the continuous group 0(2) (i.e., under rotations 
about the common center of the two circles defining the 
billiard) and is mathematically explained by the follow- 
ing result: The functions exp(±im6*), with m £ N* fixed, 
form a basis for a two-dimensional representation of 0(2). 
In order to suppress that degeneracy, it is necessary to 
break the symmetry under the continuous group 0(2). 
This can be done by folding the annular billiard along its 
diameter lying on the Ox axis. On that diameter, we can 
assume that the scalar field $ satisfies either the Dirich- 
let or the Neumann boundary condition. We then define 
two different half-annular billiards which are both desym- 
metrized versions of the annular ray-splitting billiard. By 
assuming that the modes which solve the problem defined 
by Eqs. 12151 satisfy also the Neumann (respectively the 
Dirichlct) boundary condition on the diameter, we re- 
cover the even eigenfunctions l|3*2|) (respectively the odd 
eigenf unctions (|33(l ) as well as the associated eigenvalue 
spectrum. Ea. l4Q| provides the Weyl formulas correspond- 
ing to these half-annular ray-splitting billiards. The fac- 
tor 1/2 in front of the first term of i|4U f> as well as the 
second and third terms are corrections which take into 
account the folding of the annular billiard and the bound- 
ary conditions on the fold. It is interesting to note i) the 
perimeter contribution given by ±(JVr /2ir)k that corre- 
sponds to the inner half-circle diameter which bounds the 
region of index TV by a Neumann or a Dirichlet boundary, 
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FIG. 1: The fluctuating part of the spectral counting function 
for the TM and TE theories (r/R = 2/10 and N = 3/4). 



ii) the term ±1/8 which originates from the two corners 
at the ends of the outer half-circle diameter and iii) the 
fact that the ray-splitting corners at the ends of the inner 
half-circle diameter do not provide any corrections. 



V. NUMERICAL CHECKS 

We have checked Eqs. l29l3*Tl andl4T?lfor various config- 
urations corresponding to different values of the parame- 
ters r/R, N and a by considering the oscillations around 
zero of the function AJ\T(k) — Af(k) — Af(k) which is the 
fluctuating part of the spectral counting function. All 
these numerical checks confirm the Weyl formulas ob- 
tained in Sections 3 and 4. 

In Figs, ^and |3 we present some results for the TM 
and TE theories which have been obtained for r/R = 
2/10 and for N = 3/4 and TV = 3. We have computed 
all the eigenvalues k m , n up to the frequency fc max = 120 
by solving Eq. For N = 3, the corresponding number 
of eigenvalues is around 4700 while, for N = 3/4, it is 
around 3500. 



VI. CONCLUSION AND PERSPECTIVES 

a) From our previous results obtained in the partic- 
ular case of annular ray-splitting billiards, we can infer 
rules permitting us to construct Weyl formulas for gen- 
eral ray-splitting billiards. The associated Weyl formulas 
providing the smoothed spectral counting functions are 
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FIG. 2: The fluctuating part of the spectral counting function 
for the TM and TE theories (r/R = 2/10 and N = 3). 



given in the usual form, i.e., by 



A 2 . C 



M{k) « —k 2 + — k + C, 

47T 47T 



(41) 



but now, the area term A, the perimeter term C as well 
as the constant term C must take into account the ray- 
splitting phenomenon: 

(i) A is the billiard total area weighted by the 
refraction index. For example, a piece of billiard of area 
a and index of refraction N provides a contribution to A 
given by 



N 2 a. 



(42) 



I 

(ii) £ is a sum of terms associated with the boundaries on which discontinuities in the physical properties 
occur. The contribution of a boundary of length I which separates a region of index N from a forbidden region is 
given by 



-N£ 

if we assume Dirichlet condition on that boundary and by 

N£ 



(43) 



(44) 



if we assume Neumann condition. The contribution of a boundary of length I which separates a region of index N 
from a region of index 1 is given by 



4 

TT 

47V 

7T 



E(l - N 2 ) - N - 1 
N 2 -l 



if N < 1, 



E 



N 2 



-TV- 1 



if N > 1, 



for the TM polarization, by 
"4 N 2 



7T 1 + N 2 

4 N 
n 1 + N 2 



K(l-N z ) + 



4 N 4 



TT 1 + N 2 



n(l -7V 4 ,|,1-/V 2 ) -N-l 



I if N < 1, 



A' 



N 2 



N 2 



+ 



1 



TT N(l+N 2 ) 



n 



N 4 



1 TT N 2 - 1 



N 4 



for the TE polarization and by 



4 a(l — N 2 ) 



4 a 2 N 2 - 1 



- \ / K(l-N 2 ) + - n ^-,-,1 

tt a z - 1 tt a\a l - 1) \ a z 2 



a 2 - 1 TT 



N- 



- N - 1 



N - 1 



e if TV > i, 



if N<1, 



'4 a{l~N 2 ) 
n N(a 2 - 1) 

in the general case (a ^ 1). 



K 



N 2 



N 2 



4 a{a 2 N 2 - 1) „ / 2 tt TV 2 - 1 
n 1 - a 2 ,-, 



tt 7V(a 2 -l) 



N 2 



N - 1 



if N > 1, 



(45a) 
(45b) 

(46a) 
(46b) 

(47a) 
(47b) 



(iii) The constant term C takes into account cur- 



vature and corner contributions. As far as the former is 
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concerned, it is given by 

+ wJ r W) (48) 

for a boundary curve T which separates a region of index 
N from a forbidden region, R(s) denoting the local radius 
of curvature along T. When T is a ray-splitting boundary 
which separates a region of index N from a region of 
index 1 the associated curvature contribution vanishes. 
As far as corner contributions are concerned, we simply 
note that ray-splitting corners with angle it/2 provide a 
vanishing contribution. 

We are just beginning to check the previous formulas 
for the various desymmetrized versions of ray-splitting 
sinai billiards. We obtain a very good agreement be- 
tween the theoretical formulas and the numerical data. 
This reinforces our opinion that they are exact. It would 
be very interesting to prove them rigorously but we are 
unable to do so. We have also tried to link the formu- 
las found for the TM polarization (a = 1) to the results 
that Kohler and Blumel ^(| have obtained for the scaled 
states of quantum ray-splitting billiards. We believe that 
such a link must exist but, unfortunately, we have not es- 
tablished it. 

b) In this paper, we have been exclusively concerned 



with the smooth part of the spectral counting function 
Af(k) for annular ray-splitting billiards. It seems to us 
possible to treat also the construction of the oscillating 
part of Af(k) as a canonical problem. By carefully taking 
into account Stokes phenomenon in the context of hyper- 
asymptotics [l^, |2^| , it might be possible to extract from 
g(s) all the periodic orbit contributions, even though the 
algebraic calculations involved are certainly enormous. 

c) Finally, it would be very interesting to extend our 
calculations to the three-dimensional case, having in 
mind applications to the domain of quantum optics and 
more particularly to cavity quantum electrodynamics. 
Indeed, as it is well-known, the optical properties (spon- 
taneous emission, stimulated emission, ... ) of atoms and 
molecules embedded in a cavity strongly depend on the 
density of states of the electromagnetic field. Because 
of that, Weyl formulas for cavities containing dielectric 
structures would be certainly welcome. 
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